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Abstract: We consider the operator F{u) = u' + f{t,u{t)) acting on periodic 
. real valued functions. Generically, critical points of F are infinite dimensional 

^ ! Morin-like singularities and we provide operational characterizations of the sin- 

CN I gularities of different orders. A global Lyapunov- Schmidt decomposition of F 

converts F into adapted coordinates, F{v,u) = {v,v), where w is a function of 
Q ■ average zero and both u and v are numbers. Thus, global geometric aspects of F 
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reduce to the study of a family of one- dimensional maps: we use this approach 
to obtain normal forms for several nonlinearities /. For example, we characterize 
autonomous nonlinearities giving rise to global folds and, in general, we show that 
F is a global fold if all critical points are folds. Also, /(t, x) = — x, or, more 
generally, the Cafagna-Donati nonlinearity, yield global cusps; for F interpreted 
as a map between appropriate Hilbert spaces, the requested changes of variable to 
bring F to normal form can be taken to be diffeomorphisms. A key ingredient in 
the argument is the contractibility of both the critical set and the set of non-folds 
for a generic autonomous nonlinearity. We also obtain a numerical example of a 
^ . polynomial / of degree 4 for which F contains butterflies (Morin singularities of 

^ ! order 4) — it then follows that F{u) = v has six solutions for some v. 
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Introduction 

In this paper we consider the differential equation 

u'{t) + f{t,u{t)) = g{t), (*) 

where the unknown u is a real function on §^ and the nonlinearity / : x R ^ M can 
assume a number of forms. 

Our approach is to study the global geometry of the operator 

u u' + f{t, u) 
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where the domain is either C-'^(S^) (the Banach space of periodic functions with continu- 
ous derivatives) or the Hilbert space of periodic functions with square integrable 
derivative. Ideally, we search for global changes of variables in both domain and image 
taking the operator F to a simple normal form. This goal has been achieved in previous 
occasions, starting with the seminal work of Ambrosetti and Prodi ([AP]) and its geomet- 
ric interpretation by Berger and Church ([BC]), who showed that the operator associated 
to a certain nonlinear Dirichlet problem gives rise to a global fold between infinite di- 
mensional spaces. Topological global cusps have appeared already in operators related to 
partial differential equations with a parameter ([BCT], [CDT]). Closer to the subject of 
this paper, McKean and Scovel ([McKS]) showed that the operator F for f{t,x) — (or 
more generally, for convex nonlinearities) is also a global fold, and raised the question of 
the global nature of F for /(t, x) = — x. The same question was asked by Cafagna and 
Donati ([CD]) and Church and Timourian ([CT]), who state and prove partial results for 
the more general Cafagna-Donati equation ([CD]), for which f{t, x) — ax + bx^ + cx^'''^^ 
for appropriate choices of a, b and c. In Theorem 5.1 and Corollary 5.5, we show that 
these nonlinearities indeed obtain global cusps. With some additional effort, we show that 
in Hilbert spaces the requested global changes of variables can be taken to be smooth. 

Actually, the operator F is simple enough that substantial insight into its global 
geometry can be obtained with the rather mild hypothesis of / being tame (see definition 
before Theorem 1.2). First, we construct a global Lyapunov-Schmidt type decomposition 
of F. Split a function m as a sum of a function of average zero (u) and a constant (u) and 
decompose domain and image accordingly: = ® B^. Writing the action of F as 

u = {u, u) V — {v, v), 

we show in Theorem 1.2 that, for each u, the correspondence w i— > 5 is a global diffeomor- 
phism. This provides a change of coordinates in the domain of F bringing it to (global) 
adapted coordinates 

{v, u) ^ {v, v). 

We immediately obtain that the inverse images of vertical lines under F are Gbres, curves 
foliating the domain and intersecting every horizontal plane exactly once and transversally. 
The study of F in a sense boils down to the study of its behaviour on the fibres: for example, 
f{t, x) = x^ produces a fold on every fibre and thus F is a global fold. 

Only the tameness hypothesis on the behaviour of / at infinity is necessary to obtain 
adapted coordinates: in particular, we obtain some results about the global geometry of 
F even when it is not proper. From Proposition 1.4, properness of / implies properness 
of F but, from Proposition 4.1, the converse is false. Adapted coordinates combined with 
properness make clear the possibility of definining a topological degree for F: the degree 
of F is just the degree of any of its restrictions to fibres. 

It is easy to see that (generically) 5*1, the critical set of F, is a manifold. Rather 
surprisingly, the global geometry of does not depend on the nonlinearity: generically, 
it is connected and contractible (Corollary 1.9). This follows from a more general theorem 
(Theorem 1.8 or [MST]) on the contractibility of regular level sets of a class of functionals 
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defined by integration. From contractibility, by topological arguments often using the 
infinite dimension of the spaces involved ([Ka], [Ku], [S]), there is a change of variables 
in the domain of F taking 5*1 to a closed hyperplane; in the Hilbert case, this change 
of variables can be taken to be a diffeomorphism but in the Banach case, it is merely a 
homeomorphism . 

We then proceed to study the critical points of F in detail. From adapted coordi- 
nates, ker DF has dimension 1 at critical points and imDF is then a closed subspace of 
codimension 1. This restricts considerably the possible nature of a generic critical point 
of F: it has to be an infinite dimensional Morin singularity ([M]). More precisely, after 
changes of coordinates F near a generic singularity u of can be written as 

{Z,xi,...,Xk-i,y)^{Z,xi,...,Xk-i,y'^'^^+ ^ Xiy^) 

i=l,...,k—l 

near zero, where Z is an element of an infinite dimensional space and Xi and y are real 
numbers. The integer k is the order of the singularity: folds and cusps are Morin sin- 
gularities of orders 1 and 2. Morin's classification and proof carry over to the infinite 
dimensional case by making use of a version of the Malgrange preparation theorem with 
an (infinite dimensional) parameter: this approach has been used in [CDT] to obtain a 
characterization of infinite dimensional cusps. The description of a Morin singularity is 
given more explicitly in Propositions 2.1 and 2.2 in terms of a collection of functionals 
Tiiji = 1, . . .: at a singular point of order k, the first k functionals have to be zero and 
some transversality relations have to hold. Given v, we may define a return map pv taking 
xq to xi if a (possibly non-periodic) solution u of (*) satisfies u(0) = xq, u{1) — xi. In 
Proposition 2.3, we relate the order of a singularity u to the order of contact between pF{u) 
and the identity at ■u(O). 

In the autonomous case, when the nonlinearity does not depend on t, also S^-, the set 
of critical points which are not folds, is (gcncrically) a connected contractible manifold. 
To show this, we need again Theorem 1.8 and an Lemma 3.1, stating that 5*1 and 5'2 are 
diffeomorphic to the simpler sets Si and S'2, critical and non-fold points of the simplified 
operator 

F ^ 

u^u' + J f{t,u{t))dt 

Again, contractibility yields a change of variables fiattening and ^'2. 

The functionals are rather complicated and we do not know of a simple procedure 
to decide if singularities of a given order exist for a fixed /. However, in the autonomous 
case, we describe in Lemma 3.5 a necessary (and essentially sufficient) criterion for the 
existence of singularities of order k for the simplified operator F. In the same lemma, we 
show that if F has a singularity of order k then F also does. 

In Section 4, we consider some special types of functions /: if / is either monotonic 
or convex for each value of the first coordinate t, we give a global description of the 
behaviour of F. Even though some of the results are simple or well known (from [McKS]), 
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they provide a convenient introduction to our approach of studying F fibre by fibre. Using 
Lemma 3.5, we give a criterion (Theorem 4.4) for autonomous nonhnearities to decide 
whether the operator F is a global fold. In particular (Corollary 4.5), polynomial non- 
convex nonhnearities / give rise to operators F with 82^0 but there are fast-growing 
non-convex nonhnearities for which the operator is a global fold. Also, local behaviour 
characterizes global folds (Theorem 4.6): generically, if all singularities of F are folds then 
F is a global fold. 

In Section 5, the autonomous nonhnearities / satisfy > with isolated zeros. The 
related operator F is then a global cusp: here we make full use of our techniques. In this 
case there does not seem to be an explicit description of the requested (global) changes 
of variables: their existence follows by topological arguments similar to those used in the 
study of the sets and S2- Lemma 5.4 is a global parametrized version of Whitney's 
normal form for cusps ([W]); the proof appears to be cumbersome but its main difficulty 
lies in verifying that Whitney's construction can be performed smoothly in a parameter. 
We present only a sketch of argument and we thank John Mather for helpful discussions. 

We finish the paper with an example of a different kind. The results in sections 4 and 
5 are enough to show that if /(t, x) is a polynomial in x of degree (i < 3 (with coefficients 
depending on t and non-zero coefficient of highest degree) , then the related operator F is 
a diffeomorphism, a global fold or a global cusp. In this case, thus, equation (*) has at 
most d periodic solutions. The number of solutions of (*) when / is such a polynomial was 
considered by Pugh, Lins Neto and Smale ([L]) who proved the bounds above and that 
the number of solutions may be arbitrarily large for d — A. We instead exhibit a numerical 
example of an autonomous polynomial / of degree four and a function u which is a Morin 
singularity of order four (a butterfly). This was accomplished by requesting that m be a 
root of the first four functionals Sj. By the normal form of F at a butterfly, there are 
points g near F{u) with five pre-images; one is presented. By a degree-theoretic argument, 
such a (regular) point ought to have an even number of pre-images, and we verified by 
solving the differential equation with a Runge-Kutta method that there are exactly six 
initial conditions giving rise to periodic solutions. 

1. Adapted coordinates and the critical set 

We consider the smooth nonlinear operator F : — > S° given by 

F{u){t) = u'{t)+f{t,u{t)), 

where / : x M — > M is a smooth function. Here and can be chosen in two 
different ways. In the H case, they are the Sobolev spaces B^ = = iy^(S-'^;]R) (the 
periodic absolutely continuous real valued functions with derivative in L^) and B^ = = 
L2(§1;M). In the C case, B^ = = Ci(§i;M) and 5° = C° = C°(§i;M). For notational 
convenience, inner products are to be interpreted in the sense even in other spaces. An 
interesting special situation is the autonomous case, in which / does not depend on the t 
coordinate. We denote the partial derivative of / with respect to the second variable by 

D2f. 
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Proposition 1.1 below obtains a formula for DF at arbitrary points, a description 

of the critical set SiF and a Lyapunov- Schmidt decomposition for the operator F in a 
neighbourhood of a critical point in the domain. In Theorem 1.2 we show the existence of 
a convenient global decomposition of F. 

Recall the familiar Green kernel k{x) — x — [xj — 1/2 (where [xj, following Knuth, 
is the largest integer not larger than x). If h is periodic (with period 1) then hi{t) — 
J k{s — t)h{s)ds is also periodic and h[{t) = h{t) — J h{s)ds, a function of average 0, so 
that A; is a kernel for the inverse, restricted to functions of average zero, of the derivative. 

Proposition 1.1: The derivative 

{DF{u)v){t) = v'{t) + D2f{t, u{t))v{t) 

is a Fredholm operator of index from to B^; furthermore, J D2f{s,u{s))ds is the 
unique real eigenvalue of DF{u), which is simple, with corresponding eigenvector is 

In particular, the critical set of F is 

SiF = |« e S^l J D2f{t, u{t))dt = o| . 

The subspace (l/w)-^ has codimension 1, is transversal to {w) and is also invariant under 
DF{u). Thus, the restriction 

DF{u) : gB^^ (l/w)^ C 5° 

is bijective. 

By an eigenvector of DF{u) : B^ B^ we mean a solution of DF{u)v = v' + 
D2f{tj u{t))v = Xv; by standard regularity arguments, solutions of this equation are always 
in B^. 

Proof: The formula for the derivative is straightforward. The expression for w follows 
from the explicit solution of the first order periodic linear ODE and 1/w is the only real 
eigenvector of the (adjoint) operator 

V ^ -v' + D2f{t,u{t))v. 



Let B^ ^ B^ ® (1) and = S° © (1) (the tilde denotes integral equal to 0) defining 
complementary projections 11^ and 11^. More concretely, u = H^u — {u — J u) + J u 
and u = Hgu = J u (we omit the domain of integration when it is S"^). Notice that (1) is 
always transversal to (l/w)-^ and that (w) is likewise transversal to B^. 
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A function / : S-^ x M — > R is wild at +00 (resp. —00) if 

f ds f ds 

J J max(l, supi^gi f{t, s)) ' Jj max(l, sup^^gi {-f{t, s))) 

for /= [0, +00) (resp. / = (— oo,0]); / is tame if not wild at ±cxo. Loosely, / being tame 
implies that a solution u can not go very far and come back in bounded time. 

Theorem 1.2: Let / : §^ x E — > R be a tame nonlinearity. Let F(u + u) = v + v. The 
map "if : ^ B^, \l/('u) = v + u, is a (global) diffeomorphism. 

The tameness hypothesis cannot be discarded. For f{t,u) = 27r cos(27rt) cosh^(ti), 
there are no periodic functions u such that u'{t) + f{t,u{t)) is constant. In particular, 
the point e -B" is not in the image of the map W above. Indeed, the solutions of 
u'{t) + f{t,u{t)) = are 

u = -arctanh (sin(27rt) + C) , C e (-2, 2). 

For C — 0, consider the solutions U- and u+ on disjoint domains (—1/4,1/4) and 
(1/4,3/4). Notice that U- (resp. tt+) is strictly decreasing (resp. increasing) with ab- 
solute value tending to infinity at the endpoints of the domain. 




The graph of any periodic function Ui, must cross the graphs of both u- and U-\- 
at times t- and t+, respectively, for which w'j,(t-) + fit-.u^it-)) > and 'u^(t+) + 
f{t+,Ui,{t+)) < 0. If u'^,{t) + f(t^Uy(t)) — V for all t then, from the conditions above, 
= 0. This, however, implies that Uy must equal both u- and a contradiction. 

We provide some equivalent, more geometric, readings for the rather dry statement 
of Theorem 1.2. The following diagram may be helpful: 

_ F 

U,U " V , V 




V , u 
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The map F = F o ^r-^ : S° ^ S° takes {v,u) to {v,v) = (v, (l){v,u)). Horizontal 
hyperplanes + {c} are injcctively taken by F onto sheets, i.e., hypersurfaces intersecting 
each vertical line {5} + (1) transversally and exactly once. Equivalently, the inverse images 
under F of the vertical lines {v} + (1) foliate by Gbres, i.e., curves intersecting each 
horizontal hyperplane transversally and exactly once; we denote by t„ the fibre containing 
u. Let T be the set of fibres: from transversality, we may identify T with any horizontal 
hyperplane in the domain, in particular with B^. The set of vertical lines in the image is 
naturally identified with B^ and F induces a diffeomorphism from T to B^. 

The following proof applies to both cases but certain complications are relevant only 
for the H case. 

Proof: In order to invert a vertical line and obtain a fibre, we consider the differential 
equation 

u'{t) + f{t,u{t)) = v{t) + i^, (*) 

where v e B^ is fixed and e R is a parameter. Local existence, uniqueness and continuous 
dependence on parameters hold even when v is only L^. Also, solutions cease to exist only 
by going to infinity. 

Given to and u{to) there are e > 0, and V- such that the two solutions u+ and U- 
of (*) with initial condition u{to) satisfy: 

(+) u+ either goes to +oo at some time t, to < t < to + e or satisfies u+{to + e) > u+{to). 

(— ) M_ either goes to — oo at some time t, to — e < t < to or satisfies U-{to — e) < tt-(to)- 

We discuss only (+): item (— ) is analogous. Notice that the claim is trivial in the C 
case: choose the parameter z/_|_ so that the derivative at time to of u+ is positive. Clearly, 
if (+) is satisfied by some z/_|_, it is satisfied by sufficiently positive z/_|_. 

Solve (*) for z/ = to obtain a solution uq defined on [to, to + e]. Without loss, 
uoito) > uo{to + e). Let u^^^{t) = uo{t) + e~^{uo{to) - uo{to + e)){t - to). We choose u+ 
such that the vector field (1, —f{t, u) + v{t) + i'+) always crosses the graph of Waux upwards, 
i.e., -f{t,u^^^{t))+v{t) + u+ > u[^^{t) = -f{t,uo{t)) + v{t) + e-\uoito)-uo{to + e)). 
This is clearly possible since / is continuous. □ 

Given u{0) = c, there is some z/_i_ for which the maximal solution u,y^ of (*) either goes 
to +00 at some time t, < t < 1, or satisfies u,y_^_{l) > c. This follows from the previous 
claim by a compactness argument. □ 

Given u{0) = c, there is a unique v for which (*) admits a periodic solution. Consider 
the set A'^ (resp. A~) of i/'s such that goes to +oo (resp. —oo) or satisfies u,^{l) > c 
(resp. u{l) < c). From the previous claim (and the obvious counterpart), both sets are 
non-empty. Set z/q = supA~ = inf A"*": if i/o G A~ n then Ui,^ is periodic; uniqueness 
follows from the local behaviour of the solutions. We show that fo ^ A^ implies that / is 
wild at ±cxo. We consider the case uq ^ A~ . 

If Ui^Q is defined in [0,1] with Uiyo(l) > c then continuous dependence implies that 
some open neighborhood of vo is contained in A+, a contradiction. Define tc G (0, 1] by 

lim Uuo{t) = +00. 
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If tc = 1, continuous dependence again implies that i/q is in the interior of A'^; we therefore 
have tc < I. Thus, for every M e R there exists i' < I'o and t^, > tc such that u,^{tc) > M, 
u^{ty) < c. Set 

h^{te [0, tc] I - fit, u^it)) < v{t) + v}. 

For t G we have 'U^(t) < 2'D(t) + Iv and the Lebesgue measure fi{ui,{Iy)) is bounded 
above by 2\u\ + 2\\v\\li. Define : [c, M] ^ [0, by h{s) = inf{t e [0, I u^{t) = s}. 
Even though h may have discontinuities, it is strictly increasing and then, for almost all s, 
h is differentiable with h'{s) = l/u'j^{h{s)). Let J/ = [c, M] \ Ui,{Iy): for s e J/, we have 
/i(s) ^ Ifi and /i'(s) > -l/(2/(/i(s),s)). Thus, 

/" 1 /■ (is 

tc > dhUf)) > / h'{s)ds > - / — - — --- 

c-H'y y fjj - j^^ y J -2 J^^ max(l,sup,^si(-/(^,^))) 

and therefore 



< 2tc + 2\u\ + 2\\v\\Li. 



max(l,supigsi(-/(t,s))) 
Since this estimate holds for arbitrarily large M, 

-+°° ds 



^0 



Iq max(l,supi£gi(-/(t, s))) 
A similar argument for the interval [tc, ti] yields 

ds 



< +00. 



[ 

Jo 



lo max(l,sup^^gi(/(t,s))) 
implying that / is wild at +oo. 



< +00, 



□ 



Given a, there is a unique v for which (*) admits a periodic solution u with J u = a. 
Consider all periodic solutions of (*) as curves in x R. Again, from the local behaviour 
of solutions, the curves are disjoint. By the previous claim, the union of all such curves 
contains the line {0} x R. By a similar argument applied to other lines, the union of the 
curves is x R and the curves form a continuous foliation of the cylinder by circles. Notice 
that circles in the foliation correspond to points in the fibre F~^{{v} + (1)). The integrals 
of the solutions are strictly increasing as a function of •u(O). Also, there are solutions with 
arbitrarily large (positive or negative) integrals, since the area between two curves goes to 
infinity as the initial condition of one of the curves does. □ 

At this point, we have that, given u and v, there is a unique u with F{u + u) = v + v 
(for some v). Thus, the function II^o o -P" is a bijection from any hyperplane + {a} to 
B^. This function is clearly smooth and, as discussed before the statement of this theorem, 
its derivative is always invertible. By the inverse function theorem, these bijections are 
diffeomorphisms. ■ 
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Remarks: 

1. Theorem 1.2 is the counterpart of the usual global domain decomposition found in the 
study of the equation Au = f{u), with Dirichlet boundary conditions and special resonance 
hypothesis on / (see [AP]). There, the task is simplified by the use of self-adjoint spectral 
theory. In our case, the derivative, unlike the Laplacian, is skew-symmetric, with purely 
imaginary spectrum containing 0, and the nonlinearity interacts at most with 0. 

2. The hard part in an eventual functional analytic proof of Theorem 1.2 is the properness 
(and hence, from local behaviour, bijectivity) of the function Ilfjo o F : + {v} — > H^; 
the necessary estimates seem to be simple only in the C case. 

For later use, we state as a lemma some consequences of the proof of Theorem 1.2. 

Lemma 1.3: Assume f to be a tame nonlinearity. Fibres are parametrized by average, 
i.e., the function 

is a difFeomorphism. Let Ua be the element of average a in Tuq. Then 

lim minua(t) = -t-cxo, lim maxwa(t) = — oo. 

o— >-|-oo t a— >— oo t 

Given to eS^, the function 

^ R, u{to) 

is also a diffeomorphism. 

The study of the (global and local) geometry of F thus reduces to the study of F = 
F o and therefore of : 5° — > M. The diffeomorphism is said to provide F with 
adapted coordinates, i.e., F{v,u) = (5, 0(5, tt)). This change of variables is convenient to 
the classification of critical points of F, as we shall see in the next section. Notice that we 
do not have formulae for F or \1/ and have to make do with w and $(w) = (H^g o F){u) = 
{(f) o \E')(u) = J f{t,u{t))dt (a somewhat cumbersome formula for W is given in Lemma 
1.5). 

Proposition 1.4: If f is proper then the operator F : ^ is proper. 

Notice that / : x M — M is proper if and only if \f{t, x)\ goes to infinity when {t, x) 
does. In particular, if / is proper, the restriction of F to a fibre takes infinity to infinity. 
Also, if / is proper then / is necessarily tame. 

Proof: Any compact set K C B^ is contained in the product of its (compact) projections, 
so without loss, K can be taken to be the product of a compact set K C B^ and an interval 
[— /c, k]. Since F o is of the form {v, u) ^ {v, cf){v, -u)), the compactness of the preimage 
of K under F o ^'"^ (and hence under F) follows from the boundedness of u in F~^[K) 
or the uniform boundedness of tt e F~^{K). 

In the C case, consider u at its global extrema: there, v{t) = u'{t) + f{t,u{t)) = 
f{t, u{t)) and properness of / gives us the required uniform bound. For the H case, assume 
by contradiction that there are Un G F~^{K) with Un{tn) > 2*^, where tn is the global 
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maximum of tt„, and, without loss of generality, that / is positive for large positive x. 
Consider the intervals 1^ = [tn — 1/10, If Un{t) > 2"'"^ for all t in for all sufficiently 
large n then fj^ Vn{t)dt = •u„(t^)-w„(t^-l/10)+ fj^ f{t,Un{t))dt > 1/10 min^>2"-i fit,x) 
goes to infinity with n; the norm of Vn & K is unbounded, and we are done with this 
case. Otherwise, let t'^ be the largest value in In for which 'Un(^n) = 2"'"-'^: Vn{t)dt = 

iJ'ni'tn) — '^nit'n) + ft" f{t:'^n{'t))dt > 2"^"^ for sufficiently large n and again we have a 

n 

contradiction. ■ 

Remark: There are simple a priori estimates yielding properness in the autonomous H 
case. Also, in the C case (even for non-autonomous /) easy estimates obtain Proposition 
1.4 without invoking Theorem 1.2. The analogous proof in the general H case appears 
to be considerably more elaborate and we preferred making use of the more geometric 
Theorem 1.2. 

The results above can be used to provide a simple definition of topological degree for 
the operator F in the case when / is proper: 

deg F = deg F = ^ sgn02 (w) , 
weF-i(w) 

where v is an arbitrary regular value of F. As usual, the right hand side does not depend 
on the choice of v: it is the degree of v i— > (l){v,v), a proper function from M to R. Prom 
the behaviour of at infinity (proof of Proposition 1.4), 



degF = sgn( lim /(t, x) ) — sgn ( lim f{t,x) 

Adapted coordinates give another simple characterization of the critical set: {v^u) is 
a critical point of F if and only if D2(j){v,u) = 0. Equivalently, w is a critical point of F 
if and only if D^{u)yV = where W is the tangent vector to at u given by the pull- 
back Wu = {D^{u))~^{l{^{u))) (1 is the vertical vector field consisting of the constant 
function 1 at each point). 

Lemma 1.5: Given u E and m e M, there is a unique a e E such that the equation 

to' + D2f{t,u{t))co = a, (*) 

has a (unique) periodic solution ui of average m. The function W is the only such uj of 
average 1; furthermore, W is strictly positive. 

Proof: 

The first claim follows from either solving (*) or from arguments in the proof of 
Theorem 1.2. For the element Ua of average a in t„. 



<(t) + f{t, Ua{t)) - J fit, Ua{t))dt = V. 
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Differentiating in a and setting W = -§^Ua, the equation (*) for W follows. Since J Ua = a, 
/ W = 1. The graphs of W and the constant function do not cross, implying positivity 
of W. ■ 

The following lemma introduces yet another characterization of the critical set SiF 
and, under generic hypothesis, establishes convenient transversality properties of these 
characterizations. This lemma will be essential for the more detailed study of singularities 
of F in the next section. 

Lemma 1.6: Let E^, S5, Eg : Si — > R be given by 

^a{u) = J D2f{t,u{t))dt, 

Ebiu) = J D2f{t,u{t))w^{t)dt, 

EeH = J D2f{t,u{t))Wu{t)dt. 

Then the three S 's differ by strictly positive smooth multipUcative factors. In particular, 
SiF is the zero level of each of these functionals and if is a regular value of any of the 
functionals, it is a regular value of all of them. 

Remarks: 

1. For an open dense set of functions /, is a regular value of E^. Indeed, taking deriva- 
tives as usual, is a singular value if and only if there is u e -B^ with J D2f{t, u{t))dt — 
and D2D2f{t, u{t)) = for aU teS^. 

2. In the autonomous case, is a singular value of So if and only if D2f has a double 
root. 

Proof: Here, Pi stands for a smooth strictly positive function. From the expression for 
w in Lemma 1.1, 

Thus, 

JO 



'dt. 



On the other hand. 



/o dt 



JQ 



whence 



'dt 



T,b{u) = Pi{u)T,a{u) I e Jo^ ' ^ " ^ " c 
Jo 

= P2(«)Sa(«). 
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Integrate from to 1 the differential equation describing W to obtain a = Sc(tt). 
Solving the equation, we have 

W(t) = W(0)e- /o ^^fMs))ds ^ /; D,fisMs))ds j' ^- /; D,/(r,«(r))dr^^_ 

JQ 

From = W(0), we obtain 

W(0)(1 - e-^"(-)) = Ee(«)e-^»(«) e'/o D2f{r,u[r))dr^^_ 

Jo 

Since (1 -e-^)/a; > 0, 

Sc(w) = P3(w)Sa(w), 

and we are done. ■ 

From now on, we shall always assume that / is generic in the sense that is a regular 
value of Ea; further generic properties will be required of / in Section 2 where we study 
in detail the singularities of F. 

We shall later want to use the simpler w instead of W: the following preparatory 
lemma allows for this interchange. 

Lemma 1.7: The vector fields w and W are positive multiples of each other on SiF. 
Furthermore, given u & SiF there is a neighborhood Uu ^ of u where we can write 

= aiW + Sc2;i, yV = a2W -\- TjhZ2 

for smooth real functions : f/^ — > M and smooth vector fields Zi. 
Proof: 

The first claim follows directly from the formulae for w and W when restricted to SiF 
(where = = 0). The displayed equations are consequences of the regularity of 
and Sc at S^F. ■ 

Remark: Actually, from results in [MST], in the statement can be taken to be the 
whole space B^. 

It turns out that the global geometry of is very simple, as we shall see in Corollary 
1.9. We need some preparation to state the key ingredient. Theorem 1.8. 

Let M be a smooth compact manifold equipped with a unit measure jj,. Given a 
continuous function gffc : M x M — > R'^, define Gk '■ B^ ^ to be the average of 
the related Nemytskii operator: Gk{v) = fj^ gk{m,v{m))dfx. We request that Qk admits 
continuous partial derivatives of all orders with respect to the second variable, whence Gk 
is smooth. 

Let Ilj : M'^ — > W be the projection to the first i coordinates. We say is a strong 
regular value of Gk if it is a regular value of the composition = Ilj o Gk for all z, 
l<i<k. 
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Theorem 1.8: [MST] Assume to be a strong regular value of G^. Then the levels Zi 

(1 < i < k) are contractible manifolds. Furthermore, there is a global homeomorphism S 
of taking each Zi to a closed linear subspace of codimension i; S can be chosen to be 
a diffeomorphism if = . 

The contractibility of the levels Zj essentially implies geometric triviality because of 
infinite dimension: recall that two infinite dimensional separable Hilbert manifolds are 
diffeomorphic if their homotopy groups coincide ([Ku]) and that all infinite dimensional 
separable Banach spaces are homeomorphic ([Ka]). 

In the next corollary we have k = 1; Theorem 1.8 in its generality will be convenient 
in Section 3. 

Corollary 1.9: Assume that is a regular value of E^. Then Si is connected and 
contractible. Furthermore, there is a global homeomorphism S of B^ taking Si to a 
closed linear subspace of B^ of codimension 1; S can be chosen to be a diffeomorphism if 
B^ = H\ 

2. Morin theory 

Morin classified generic singularities of functions from to whose derivative has 
kernel of dimension 1 ([M]). The first step in Morin's proof makes use of the implicit 
function theorem to write such a singularity at the origin in adapted coordinates, i.e., in 
the form 

(x,y) ^ (x,//(x,y)), x= {xi, . . . ,Xn-i) G y eR, 

after composing with suitable diffeomorphisms in the neighborhoods of zero in both domain 
and image. Morin's central result is that such singularities are classified by their order, a 
Morin singularity of order /c is a point (x, y) for which 

(a) D2n{x, y)^ ■■■ = //(x, y) = 0, 

(b) D2'+V(x,y)7^0, 

(c) the Jacobian D{D2iJi, . . . , -D2^^/x)(x, y) is surjective. 

Set Sk = {(x, y)|L)2A*(x, y) — ■■■ — D2iJi{^,y) = 0}. Thus, Sq is the domain. Si is 
the critical set {(x, y)\D2n{x, y) — 0} (consistently with previous notation). Also, a Morin 
singularity of order i belongs to 5"^ if and only if i > k. In a neighborhood of a Morin 
singularity, the sets Sk stratify the domain: the sets are nested and Si is a submanifold of 
codimension i. Notice that a point (x, y) E Sk — Sk+i is a Morin singularity (of order k) 
only if condition (c) above holds. 

Composing by appropriate diff'eomorphisms in the domain and image, a Morin singu- 
larity in dimension n and order k acquires the normal form 

(a;i, . . . , Xn-i, y) (xi, . . . , x^-i, 1/''+^ + xiy''~^ H h Xk-iy). 

Morin singularities of order 1, 2, 3 and 4 are called, respectively, folds, cusps, swallowtails 
and butterfiies. 
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We shall need an equivalent classification for singularities of functions between infinite- 
dimensional spaces. Let G : Zi — > Z2 be a smooth map between Banach spaces so that 
DG{zq) is Fredholm operator of index and kernel of dimension 1. Again, after changes 
of variables in the domain and image, we may assume G near zq to be written in adapted 
coordinates as 

G:XxR^XxR, 

(x,y) ^ (x,//(x,y)) 

and if conditions (a), (b) and (c) above hold, the same normal form applies for an appropri- 
ate splitting X = M'^"^ ® X' — we then call zq a Morin singularity of order k. The proof of 
this last fact follows Morin's ([M]), making use of a parameterized version of Malgrange's 
preparation theorem, the parameter taking values in a Banach space (see [CDT]). 

We already saw in the previous section that the composition F = Fo^f~^isin adapted 
coordinates: 

F : SO© (1) ^ SO© (1). 

{v,u) >—>■ (^v,v = (f){v,u)) 

From the previous paragraph, a point is a Morin singularity of F (or F) of order k if and 
only if conditions (a), (b) and (c) hold for fi replaced by (p. This criterion, however, can not 
be used directly since we have no formula for (p; we rephrase it in terms of $ = 0o^' and w. 
Following the usual notation, we write w$, for the Lie derivative D^{u) ■ Wu- The following 
result shows that we may substitute w for W (alternative generators for the kernel of DF 
over Si) — a fact which in finite dimension would be unsurprising. 

Proposition 2.1: The point u E is a Morin singularity of order k for F if and only 
if: 

(a) w^{u) = ■■■ = w''^{u) = 0, 

(b) w''+^^{u) ^ 0, 

(c) D{w^, . . . , w'^~^^){u) is surjective. 

Proof: Consider in 5°© (1) the constant vertical vector field 1, consisting of the constant 
function 1 at each point. In the notation we just introduced, = 1^. In terms of the 
pull-back W{u) = (DW(w))~^(l(W(w))), the conditions for w to be a Morin singularity of 
order k are: 

(a') W^{u) = ■■■ = yV''^{u) = 0, 
(b') W'=+i$(w) 

(c') D{W^, . . . , >V^-i$)(m) is surjective. 

We are left with showing that we can substitute W by in these conditions. Notice 
first that w$ = and >V$ = S^,, proving the case k — (regular points). From now 
on, we assume u e SiF and write, making use of Lemma 1.7, w = aiW + {W^)zi and 
W = a2W + {w^)z2 in a small neighborhood Uu of u. 
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For each k, the ideals in C°°{Uu, IR) generated by w^, . . .w^<^ and >V^», . . .W^^ are 
equal. Assuming by induction that the result holds for A; — 1, 

= (aiW + + ■ ■ • + bk-iW''-^^) 

= ai(W6i)(W$) + aibiW^^ + ■ ■ ■ ai{Wbk-i){W^~^^) + ai6fe_i>V^$+ 
(>V$)2i(6i>V$ + bk-iW^~^^) 

which is clearly in the ideal with generators >V^», . . . W'^^*, proving one inclusion; the 
opposite inclusion is analogous. 

The equality of the two ideals with k generators implies 

w''^ = bkW''^ + ■■■ + 61 W$, 
W^^^ = Ckw'^^ H h ciw^, 

for smooth functions bi and Cj where bk and Cfc are non-zero. The equivalence between the 
conditions (a) and (a') or (b) and (b') is clear. The third equivalence follows from repeated 
use of the simple fact that the spans of D{gi{u)g2{u)) and D{gi{u)g2{u) + a{u)gi{u)) 
coincide for points u such that gi{u) — g2{u) = {a being a smooth real function). ■ 

Proposition 2.2: For 

Si(w) = J D2f{t,u{t))dt, 
S2(w) = J Dy{t,u{t))w{t)dt, 
^3{u) = J Dlf{t,u{t))w\t)dt, 

E4(«) = J Dif{t,uit))w\t) -2Dlf{t,uit))w\t) (^j^ Dlf{s,u{s))w{s)d^dt, 
S5(«) = j Dlf{t,u{t))w\t)-5Dif{t,u{t))w^{t) (^^ Dlf{s,u{s))w{s)ds^ + 

bDlf{t,u{t))w\t)(^j^ Dlf{s,u{s))w{s)^ dt, 
we have, for k = 1, . . . , 5, 

Sk = {ue = 0,i = 1, .. 

Furthermore, for k = 1, . . . ,4, u is a Morin singularity of order k if and only if 

(a) Ej(tt) = 0, i = 1, . . . ,k, 

(b) Sfc+i(«)7^0, 
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(c) the derivative DIl{u) of the function 

u . . . , Sfe_i(tt)) 

is surjective. 

Proof: From Proposition 2.1, we want to compute w'^^{u). The expressions for T,k{u) 
follow from repeated integration by parts, discarding elements in the ideal generated by 
w^{u), . . . ,w''~^^{u) and non-zero multiplicative factors. ■ 

In particular, w is a fold point if and only if Ei(ti) = and S2(ti) ^ 0. Also, u is 
a cusp point if and only if 'Ei{u) = E2(tt) = 0, Tis{u) ^ and DYi\{u) ^ 0. Clearly, 
-DEi(u) ■ V — J D2f{t,u{t))v{t)dt, and hence -DEi(u) = if and only if the function 
D^fit, u{t)) is identically 0. 

There is a simple relationship between Morin singularities and the return map. Given 
/ : §^ X R ^ R and V e the return map p„ : / ^ R, / C R, sends u{0) to u{l) if 
tt : [0, 1] — > R satisfies 

u\t)+f{t,u{t)) = v{t). 

Here / is the maximal domain, i.e., the set of initial conditions such that the solution u 

extends to t = 1. 

Proposition 2.3: For u e and v = F{u), pv{u{0)) = u{0). Also, u E Si if and only 
if p;(w(0)) = 1 andueSk if and only if p'^{u{0)) = 1 and pi"\u{0)) ^ for i = 2, . . . ,k. 

Proof: Consider the fibre Tu through u. Let Ua be the element of average a of r„ and 
u = Uao- by Lemma 1.3, is smoothly parametrized by a. Let g{a) — ^{ua), the average 
of F{ua); clearly, g^'^\a) = >V'^$(ua) and therefore u e S'^ if and only if g''''\ao) = 
for i = 1, . . . ,k. Let u'^ be the element of with tt'^(O) = c: by Lemma 1.3, is also 
parametrized by c. Let h{c) — ^{u'^) —v = F{uc) — F{u): by the chain rule, u E if and 
only if (w(0)) = for i = 1, . . . , fc. 

For c, 6 e R, let /3(c, b) = py+b{c)'- that is, if U (c, b, t) satisfies 

DsU{c,b,t) + f{t,U{c,b,t)) = v{t) + b, C/(c,6,0) = c (*) 

we have /9(c, b) = U{c,b,l). The periodicity of it^ yields /3(c, h{c)) = c. Points in the curve 
(c, h{c)) thus correspond to points in and the largest k for which u E Sk is the order of 
contact between this curve and the horizontal axis (c, 0) at the common point (14(0), 0). 

Differentiating (*), 

D3D2U{u{0), 0, t) + D2f{t, u{t))D2U{u{0),0, t) = 1, £>2C^(«(0), 0, 0) = 

and, by explicitly solving for D2U, we obtain D2U{u{0),0, 1) = i3>2/3('u(0), 0) > 0. Thus, 
G{c,b) — {c,(3{c,b)) is a local diffeomorphism near (u(0),0) taking the curve (c, /i(c)) to 
the diagonal (c, c) and the horizontal axis to (c, p„(c)). The order of contact between the 
curves is preserved by G and we are done. ■ 
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3. The autonomous case 

This section is dedicated to a number of special properties of the autonomous case, 
when / depends on x only. Notice that this implies that / is tame. 

The sets Sk are described by the rather complicated formulae S/-; in the autonomous 
case it is convenient to consider the simpliiied operator 



F : ^ 

u ^ u + J f{u{t))dt 



whose critical strata Sk are far easier to handle but still convey significant information 
about F and Sk- 

Since F is already given in adapted coordinates, straightforward application of Morin's 
characterization obtains 



Sk^{veB^\ J %{v{t)) = 0}, 



where ■jkix) is the A;-dimensional vector (/'(x), . . . , f^'^\x)) Notice that, from Lemma 1.6, 
= 5*1. However, §2 is usually different from ^'2: the following lemma relates both sets 
and is a key ingredient in Section 5. 

Lemma 3.1: Let / : M ^ M be a smooth function, F : B^ B^ be the operator 
{F{u)){t) — u'{t) + f{u{t)). Then there exists a global diffeomorphism of B^ taking Si to 
itself and S2 to 82- 

Proof: The diffeomorphism has the form u ^-^ v = u o a, with inverse v ^-^ u = v o 
where a and /? are orientation preserving diffeomorphisms of S-*^ fixing 0: such 

compositions take functions in B^ to functions in B^. Clearly, /? = a~^. We need the 
following characterizations, which follow easily from Lemma 1.6 and Proposition 2.2: 

5i = |w e J f'{u{t))w{t)dt = o| , 

ueB^ J f\u{t))w{t)dt = J f"{u{t))w{t)dt = 

We first obtain /? from u: 

, , C w(a)d(T 
lo w{a)da 

is clearly a diffeomorphism {w > 0). For any continuous function g, a change of 
variables gives 

/ g{u{s))w{s)ds ^ ^ [ g{v{t))dt^O, (*) 

Jo Jo 

where u = v o (3. Thus, if tt e 5 then v E S and if u E C then v E S2- The smooth 
dependence of /? on tt is obvious. 
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To show invertibility of the map m i— > we obtain a from v. \i a satisfies 

oiit) = , a(0) = 

A-ljlnv{T))dT 

for an arbitrary positive constant standard algebra shows that 

a'{t) = J exp f'{viT))a'iT)dT^ , a(0) = 0. (f) 

This again imphes the equivalence (*) (where, of course, v = uoa) and hence v E S (resp., 
S2) implies u E S (resp., C) provided a(l) = 1. We have to show that for each v there is 
a unique A with a(l) = 1 and that the dependence of A on v is smooth. 

Let h{t) = f'{v{T))dT; h is with h{Q) = 0. The function a is defined in [0, 1] if 
A > maxf h{t). From 

= i A^y 

the derivative of a{l) with respect to A is strictly negative. When A tends to infinity, 
a{l) becomes small and when A approaches max^ h{t), a{l) tends to infinity. This settles 
existence and uniqueness of the required A; smoothness follows from the implicit function 
theorem applied to the smooth function {v,A) 1— > a{l) and the fact that the derivative 
with respect to A is not zero. 

It remains only to show that the two smooth maps constructed above are the inverse 
of each other. Consider the sequence of maps 

a _l /3 ~ n-1 

V — m — V o a — >v — uo [j 
We have, for positive constants Ci and C2, 
13' {s) = Ciw{s) 



^Ci exp ^- ^ ^ f'{v{T))a'{T)dT^ (from r = a'^a)) 

(t) C2 
a'(a-i(s)) 

= C2{a-y{s) 

and, since a{0) = /?(0) = and a{l) = /?(1) = 1, it follows that = (3 and v = v. 
Similarly, for 

/3 n—\ a ~ _i 

u — w = uo [3 — m = V oa , 
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we have 



a'it) = Csexp(^l^ f'{v{T))a'{T)di 

^=C('^)c3exp I r^'^ f'{v{a-\a)))da 



'0 



w{a{t)) 
(3'{a{t)y 

hence {(3 o a)' is a constant and the result follows. ■ 

As with S\ — Si, the global geometry of the sets Sk is very simple, as we learn from the 
following application of Theorem 1.8. We say that the simplified operator F is k-regular 
if is a strong regular value of 

E : ^ 

u i-> (Si, . . . , Sfc) 

Corollary 3.2: Assume F is k-regular. Then there is a global homeomorphism S of 
taking each Si (1 < i < k) to a closed linear subspace of B^ of codimension i; S can be 
chosen to be a diffeomorphism if B^ = . 

Combining Lemma 3.1 and Corollary 3.2 we have: 

Corollary 3.3: Assume F is 2-regular. Then there is a global homeomorphism S of 
B^ taking each Si, i = 1,2, to a closed linear subspace of B^ of codimension i; S can be 
chosen to be a diffeomorphism if B^ = . 



It seems hard to give an operational criterion to decide for larger k even whether S^ 
is non-empty. We now present a partial criterion. 

Definition 3.4: Let / : M ^ M be a smooth function; f is said to be /c-good if 7^ never 
vanishes and the image of any open interval by 7/5 is not contained in a hyperplane through 
the origin in M^. 

Generic smooth functions are /c-good, as well as generic polynomials of fixed degree at 
least k. It is easy to see that if / is fc-l- 1-good, then the simplified operator F is fc-regular. 

Lemma 3.5: Let / : M — > M be a k-good function and F be the related operator. Then 
Sk ^ if and only if e IR'^ is in the interior of the convex hull of the image of 7fc. Also, 
Sk ^ implies Sk 7^ 0- 

In section 6, we give an example of a polynomial of degree 4 (and hence for which 
S4 = 0) having ^4 7^ 0. 
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Proof: 

Clearly, if 0, is in the convex hull of the curve 7^ (but not on the curve itself). 
If is on the boundary of the convex hull, by a standard support theorem ([G], pg. 12), 
the image of is contained in a closed half-space defined by u{p) > for some linear 
functional u. Since / is /c-good, for any non-constant v G H^, u{J 7A;(f)) > and Sk = 0. 

Conversely, assume in the interior of the convex hull of the image of 7/s. By Steinitz's 
theorem ([G]), there are points 7fc(a;j), j = 0, . . . , 2A; — 1 such that is in the interior of 
their convex hull. For e e (0, 1) and aj > 0, j = 0, . . . , 2/s — 1, ^ = 1 — e, consider a 
smooth function v^^a of period 1, defined as follows. We split the domain [0, 1] into intervals 
/q, Jo, ... , l2k-ij J2k-i of lengths uq, e/2/c, . . . , 0.2^-1, e/2/c; inside /j, v^^a is constant equal 
to Xj and inside each Jj, v^^a is the appropriate affine transformation of a fixed smooth arc 
joining two steps. As e tends to 0, Ve,a approaches a step function. Let (/)(e, a) = J ^k{ve,a)' 
this function is affine in e and a (i.e., linear plus constant) and to show that Sk is non- 
empty, we need to find a zero of 0. The function extends continuously to e = and is 
then in the interior of the image of the simplex spanned by the Oj's: there exists therefore 
a straight segment parametrized by small positive values of e on which (p is zero and Sk 
is thus non-empty. More, for a fixed small eg, take a /c-subspace V of M^'^ such that, for 
a & V, 0(eo, a) is surjective. The image under of a small sphere around the origin in V 
is some ellipsoid containing the origin in its interior. 

Now, assume S^ 7^ 0. Use the space V and the function to obtain r > and a 
function H ■.M'' ^ with 

J %{H{s){t))dt = rs 

for s e B'^ where B'^ C M'^ is the unit ball. Define the N -replicator to be the isomorphism 
Rn : B' ^ Rn{B') C B\ {RN{u)){t) = u{Nt), i = 0, 1. Clearly, (i?iv(w))' = NRn{u'). 
We claim that given e > there exists N such that 

\{w^Rn{H{s))),...,w''^Rn{H{s)))) -rs\<e 

for all s eM^ and the proof is completed by a standard degree theory argument. 

At this point it is convenient to make explicit the dependence of w = w{u) in terms 
of u. From Proposition 1.1, 

{w{RN{u))y{t) + f'{RN{u{t)))w{RN{u)){t) = Xw{RN{u)){t), A = y f'{RN{u){t))dt 

and therefore A is the same for u and Rn{u). Define wn{u) by Rn{wn{u)) — w{Rn{u)) 
so that 

and, from the formula for w in Proposition 1.1, WN{u){t) = {w{u){t)y^/^\ Since $(w) = 
/ f{u{t))dt we have ^{Rn{u)) = $(u) and 

{w^{Rn{u)), . . . , w''^{Rn{u))) = {wN^{u), . . . , w%^{u)) . 
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The sequence {wn) of vector fields tends to the constant vector field 1 (i.e., the constant 
function 1 at every point u) in the C^-metric (for any n). Also, 

{1^Rn{u)),...,1''^Rn{u))) = J %{u{t))dt, 

proving the claim. ■ 

Remarks: For k = 2, A;-goodness may be weakened to /' and /" having no common zeros. 

It also follows from the above argument that D{{w^, . . . ^w^^) o Rj^ o H) tends to 
the identity matrix when goes to infinity, establishing condition (c) in Proposition 2.1. 
Condition (b) follows from the additional hypothesis of / being (A; + l)-good, thus proving 
the existence of Morin singularities of order k. 



4. Some examples 

In this section, we describe the global geometry of F : ^ for several special 
classes of smooth functions / : §^ x M — > M. The first rather simple proposition illustrates 
the use of fibres and adapted coordinates in three technically different scenarios. 

Proposition 4.1: 

(a) If f is proper and D2f{t, x) > then F is a diffeomorphism. 

(b) Assume D2f{t,x) > 0. Let f±{t) = lim.x->±oo f{ti x). Then F is a diffeomorphism 
from B^ to the horizontal strip 

S^{veB^\ J f-{t)dt < J v{t)dt < J f+{t)dt}. 



(c) If f is proper and strictly increasing in the second variable then F is a homeomorphism. 
Proof: Notice that the hypotesis on each item imply that the nonlinearity / is tame. 

(a) Recall that, by Theorem 1.2, B^ is foliated by fibres, B^ is foliated by vertical lines and 

there is a diffeomorphism between the space of fibres and the space of vertical lines. From 
the characterization of Si in Proposition 1.1, we see that F has no critical points. Thus, 
F takes each fibre strictly monotonically to its related vertical line in B^ — by Proposition 
1.4, F is actually a diffeomorphism from fibre to vertical line and the result follows. 

(b) As in the previous item, fibres are bijectively taken to open subintervals of vertical 
lines. More explicitly, a function u in the fibre (H^o o F)-^(v) is taken to u' + f(t, u(t)) = 
V + J f{t,u{t))dt and the extremes of the image of this fibre are lima^_|_oo / f{t,Ua{t))dt 
and lima^-oo / f{t,'iJ'a{t))dt, where Ua is the element of average a in the fibre. From 
Proposition 1.3, 

lim miii Ua{t) = +oo, 

o— >+oo t 

lim maxtt(j(t) = — oo. 

a— > — oo t 
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Thus, 

lim / f{t,Ua{t))dt= I f+{t)dt, 
^Hrn^y" f{t,Ua{t))dt = J f-{t)dt, 

and the result foUows. 

(c) By properness (Proposition 1.4), it suffices to prove that F is strictly increasing on 
each fibre (notice that F restricted to fibres may have critical points where WF — 0). Let 
uq < ui be elements of the fibre (H^o o F)~^{v) so that 

F{ui){t) = u[{t) + /(t, Ui{t)) = v + Vi. 

Integrating in t e and using the monotonicity of /, we obtain vq < wi, concluding the 
proof. ■ 

Remarks: 

1. A more standard proof of (a), without making use of Theorem 1.2 (and the consequent 
fibre-sheet-adapted coordinates vocabulary), could be as follows. As before. Si ~ and 
from Proposition 1.4, F is proper. Since is simply connected, by covering space theory, 
F is a diffeomorphism. Notice that this argument does not extend easily to the other 
items. 

2. Similar results and proofs hold if instead D2/ < 0. 

3. In item (b), if both J f^{t)dt and / f-{t)dt diverge, then F is a global diffeomorphism. 
In particular, F may be proper even if / is not. 

Theorem 4.2: If / is proper and -Dl/l*) a;) > then F is a global fold. 

We call an operator G : ^ a global fold if there exist diffeomorphisms Si : 
^ R X 5° and So : 5° ^ M X S° such that (Sq o G o E^^){x,v) = {x^,v), for all 
(x^v) E MxS°. Similarly, we call G a topoJogicai global fold if there exist homeomorphisms 

Si as above. 

Proof: From D|/(t,a;) > 0, we conclude that D2{t,x) is strictly increasing in x for any 
fixed t and hence that Sa(w) = J D2{t,u{t))dt is strictly increasing on fibres. Thus, each 
fibre contains a unique critical point uq and, for arbitrary u- and in the same fibre as 
Uq satisfying tt_ < < tt+, we have Ea(tt_) < and T,aiu+) > and, from Lemma 1.6, 
W$(w_) = Sc(w-) < and W^{u+) = Sc(w+) > and the restriction of $ to a fibre is 
a global fold from M to M. Thus, on each fibre, we have diffeomorphisms S^ as above and 
the problem is whether such diffeomorphisms can be chosen so as to depend smoothly on 
the fibre. 

In adapted coordinates, wc must define ^1 : x M — > M and ^0 : -^^ x 1^ — >^ so that 
the vertical columns of the diagram below are diffeomorphisms and the diagram commutes, 
i.e., 

^0(5,0(5,^)) = (^i{v,u)f. 
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We construct the explicitly. For each v, let ay be the unique critical point ofu\-^ (/»({;, u). 
Clearly, ay and its image bj, = (f>{v, ay) depend smoothly on v. Also, write (p{v, u) — by = 
{u — ayY g{v^u). From the previous paragraph, (7 is a smooth positive function. Set 
io{v, v) = V - by and ^i{v,u) = (u - ay)y/g{v,u). ■ 

Remarks: 

1. McKean and Scovel ([McKS]) studied this scenario with a different set of fibres for 
and its image under F. Our choice of fibering by vertical lines and inverting them 

under F to get a fibration oi B^ is more helpful in our examples. 

2. Slight variations (as in Proposition 4.1) are possible and can be handled similarly; we 
omit the tedious details. 

In the autonomous case, theorem 4.2 admits a partial converse: 

Theorem 4.3: Let f : M. ^ M. be a smooth function which is both 2- and 3-good, with 
lima;^±oo f{x) — +00. If is in the interior of the convex hull of the image of 72 then F 
has cusps. Otherwise, F is a (differentiable) global fold. 

Recall that 72 (t) = {f'{t),f"{t)). Notice that once F is known to have a cusp, from 
the normal form we have image points with three regular pre-images near the cusp; since 
deg(F) = 0, such points have at least one additional pre-image. 

Proof: If is in the interior of the convex hull of the image of 72, 5'2 7^ by Lemma 3.5. 
We must prove that some points u in 5*2 are differentiable cusps, i.e., satisfy Ti^^u) 7^ and 
DEi(ti) ^ 0. As remarked after Proposition 2.2, D'Ei{u) = only when f"{u{t)) = for 
all t, which implies, given 2-goodness, that u is constant equal to a root of /". Again from 
2-goodness, /' and /" have no common roots and T,i{u) = J f'{u(t))dt ^ 0, thus u is not 
in the critical set. It remains to verify that we may choose u so that S3 (it) ^ 0. From 3- 
goodness, the curve 73(t) = {f'{t),f"{t),f"'{t)) docs not intersect the origin. The convex 
hull of the image of 73 meets the vertical axis and must contain points (0, 0, A) distinct 
from the origin, otherwise the image of 73 would have to be contained in a hyperplane, 
contradicting 3-goodness. Imitating the proof of Lemma 3.5, we may construct u with 
Si(tt) = T,2{u) = 0, T,^{u) A, which is the required cusp. 

Conversely, assume not to be in the interior of the convex hull of the image of 72. 
Again from Lemma 3.5, all critical points are folds. Notice that if m is the minimum of /, 
f'{ni) = and, by 2-goodness, f"{m) > 0. Thus, the intersection of the convex hull with 
the second axis consists of points with non-negative second coordinate. From the proof of 
Lemma 3.5, Si(u) = now implies Tj2{u) > 0; thus, the fold points have concavity upwards 
in the restriction of F to each fibre. From the behaviour of / at infinity, F has exactly one 
fold point per fibre. The rest of the argument is similar to the proof of Theorem 4.2. ■ 

Corollary 4.4: 

(a) Let / : M — > M be a 2- and 3-good polynomial of even degree and positive leading 
coefRcient. If f" assumes both signs then the operator F has a cusp and there are points 
in the image of F with four pre-images. 

(b) There are non-convex smooth proper functions / : R — > M for which F is a global fold. 
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Proof: To prove (a), notice that, by hypothesis, there is a point in the lower half-plane 

(/" < 0). Also, since / is a polynomial, the argument of 72 (t) approaches from above 
(resp., TT from below) when t goes to +00 (resp., —00). This suffices to show that is in 
the interior of the convex hull of the image of 72 (t)- 

As for (b), we deviate from the previous argument by considering functions for which 
/' and /" are comparable for large x, such as cosh a;. More precisely, let / be such that 
/" coincides with cosh t outside a small interval of large positive numbers in which we 
subtract from cosht a narrow positive bump — if the bump is sufficiently high and narrow, 
/" changes sign, f'{x) > for all positive x and is not in the convex hull of the image of 
72- ■ 

We can also characterize global folds in a purely local way: 

Theorem 4.5: If f : x M. ^ R is proper, is a regular value of Si and all singularities 
of F are folds, then F is a global fold. 

Proof: For the operators F being considered, there is a sign associated to each fold: it 
is the sign of S2, which can also be interpreted as saying whether the concavity of the 
restriction of F to fibres points up or down. This splits the set of folds into two open 
subsets. Since by hypothesis S2 — and by Corollary 1.9 S'l is connected it follows that 
one of these sets is empty. In other words, all folds are concave upwards, say, and the 
restriction of F to any fibre thus has at most one critical point. The result follows by 
properness and juxtaposition of fibres as in Theorems 4.2 and 4.3. ■ 

5. Global cusps 

The results in the previous section describe the global geometry of F when / satisfies 
D^f > or D^f > 0. In the autonomous case, we are able to handle another kind of 
nonlinearity. 

Theorem 5.1: Let / : R — > R be a proper function such that 

(a) fix) > 0, 

(b) f"'{x) has isolated roots, 

(c) f'{x) assumes both signs. 

Then F is a topological global cusp; in the H case, F is a smooth global cusp. 
A similar result holds if (a) is replaced by f"'{x) < 0. 

The scenarios of the previous section arc simple enough to allow for rather explicit 
global changes of variable to normal form. This is partly due to the fact that restrictions 
of F to arbitrary fibres have similar behaviours: in a sense, we may consider fibres in- 
dividually. For the operators in the statement of Theorem 5.1, instead, such restrictions 
vary according to the fibre, as illustrated in Figure 5.1: we must therefore treat them 
collectively. In the process, explicitness is lost: from the theorem (and its proof), in the 
H case the domain and image are foliated by smooth surfaces, diffeomorphic to M^, which 
are in turn foliated by fibres. More, F takes surfaces to surfaces and, on each surface, F 
is a global cusp. Still, we have no idea how to exhibit such foliations: they are shown to 
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exist by topological methods in Hilbert manifolds ([Ku], [MST]). For Banach spaces, we 
use an additional existential argument depending on the homeomorphism of all infinite 
dimensional separable Banach spaces ([Ka]). 

An operator G : is a smooth global cusp if there exist diffeomorphisms 

Si : ^ HxR'^ and So : H xR"^ such that {EooGoE~^){Z,x,y) = {Z,x,j/ + xy), 

for all {Z,x,y) e H xM.'^, where H is a, separable infinite dimensional Hilbert space. 
Similarly, we call G : — > C° a topological global cusp if there exist homeomorphisms 
Si : ^ X R2 as above. 

The following lemma is an exercise in integration by parts. 

Lemma 5.2: Let ^ : R — > M be a smooth function. Given a < b, we have 



An additional topological ingredient is an infinite-dimensional version of the corollary 
to Theorem 1 in [S]. The proof is similar for Hilbert spaces and makes use of results in 
[BH] for Banach spaces (see [MST] for additional information). 

Lemma 5.3: 

(a) Given a contractible connected smooth submanifold H' of codimension 1 of a separable 
Hilbert space H of infinite dimension, there is a diffeomorphism of H to itself taking H' 
to a closed subspace of codimension 1. 

(b) Let B' be a closed subset of a separable infinite-dimensional Banach space B. Assume 
that B' is connected, contractible, and a bicoUared topological submanifold of codimension 
1 in B. Then there is a homeomorphism from B to itself taking B' to a closed subspace 
of codimension 1. 

Actually, a contractible connected closed bicoUared topological submanifold of codi- 
mension 1 always splits the ambient space in (exactly) two components. 

Finally, we make use of a canonical construction to bring planar global cusps to normal 
form. A sketch of proof is given at the end of this section. 

Lemma 5.4: Let Z be a topological space of parameters. Let G : Z xM? ^ Z xM? be 
a continuous function of the form 





G{Z,x,y) = {Z,x,gz{x,y)) 



with the following properties: 
o For aU Z eZ, gz (0,0) = 0. 

o For all Z & Z and a; e R, limy^±^ gz{x, y) = ±oo. 

o For all Z e Z and x <0 the function y i— > gz{x, y) is strictly increasing. 
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o There exist continuous functions m, M : Zx[0, +00) R with m{Z, 0) = M{Z, 0) = 
such that, for all Z & Z and x > the function y gz{x, y) is strictly increasing in 
{—00, m{Z, x)] and [M{Z, x), +00) but strictly decreasing in [m{Z, x), M{Z, x)]. 

(a) There exist homeomorphisms Wd and Wi of Z x keeping the Z and x coordinates 
fixed such that 

{WioGo W^'){Z, X, y) = (Z, X, y^ + xy). 

(b) If Z is a smooth Hilbert manifold and G is smooth with 

D2gz{0, 0) = Dlgz{0, 0) = 0, Dlgz{0, 0) > 0, DiD2gz{0, 0) < 

then Wd and Wi can be taken to be diffeomorphisms. 

Proof of Theorem 5.1: 

We first classify tfie singularities, next we study the behaviour of the restriction of F 
to fibres and finally obtain global results by juxtaposing fibres. 

All critical points of F are folds or cusps. Clearly, E3 is always positive and we only 
have to check transversality when Ei = E2 = 0, i.e., we have to verify that DEi at 
such points. Since Si = J f'{u{t))dt, 



DJ:i{u) J f"{u{t))v{t)dt; 



we show that the function f"{u{t)) is not identically zero. First, there is a unique xq for 
which f"{xo) = 0; from (c), we have f'{xo) < 0. Thus, the function f"{u{t)) is identically 
zero only if u{t) = xq for all t but then J f'{xt)dt < and u ^ Si. □ 

We now consider F restricted to fibres. From Theorem 1.4, such restrictions take ±00 
to ±00. Also, regular and fold points of F are regular or fold points of the restriction. 
Furthermore, the restrictions are locally increasing at cusp points. Indeed, S3 > and, at 
cusp points, S3 is a positive multiple of >V"^$, the third derivative of the restriction. 

The operator F has at most two critical points per hbre. Let ui < U2 be two critical 
points of F in the same fibre: 

J f'{ui{t))dt = J f'{u2{t))dt = (*) 

and 

u[{t) + f{m{t)) = u',{t) + f{u2{t)) + C (t) 
for some constant C . We first show that C > 0. From Lemma 5.2, 

/-(..w)+/'fa(o)- ^'";';»"y >o 

U2{t) - Ui{t) 
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and, integrating and making use of (*), 



/ 



fiu2{t))-f{u,{t)) 
U2{t) - Ui{t) 



dt < 0. 



From (t), 



{U2-Ui){t) U2{t)-Ui{t) 



+ 



U2{t) -Ui{t) 



C 



= 



and, integrating, we obtain 




dt 



> 



•U2(t) - Ui{t) 



implying C > 0. This means that the restriction of F to fibres, if further restricted to 
critical points, is decreasing: ui < U2 implies F{ui) > F{u2)- It follows that, if there 
are at least three critical points on a fibre, F is decreasing near the second one: from the 
previous paragraph, this second critical point can thus be neither a fold nor a cusp. □ 

At this point we already know that regular values have one or three pre-images, images 
of folds have two pre-images and images of cusps have a single pre-image (as in [CD] and 
[CT]). But we have more: F restricted to a fibre is topologically equivalent to one of the 
three graphs on Figure 5.1. Fibres containing a cusp, on which F behaves as depicted in 
(b), split the space of fibres into two subsets, on which restrictions of F behave as in either 
(a) (no critical points) or (c) (two folds). Let T = be, as before, the space of fibres: 
we have the natural partition !F = !Fa^ into sets of fibres of types (a), (b) and (c), 

respectively. 



The sets Ta and Tc are open in T; Th is closed. Let uq G J^c'i there exist elements 
ui < tt2 in the tto-fibre with F{ui) > F{u2). The hyperplanes parallel io passing 




(a) 



(b) 



(c) 



Figure 5.1 
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through ui and U2 transversally intersect fibres. In particular, fibres sufficiently close to 
the MQ-fibrc contain point u[ < u'2 in these hyperplanes for which F{u']) > F{u2), proving 
the openness of J-'c- 

Assuming by contradiction that JF„ is not open, let Un be a sequence of critical points 
whose corresponding fibres converge to the fibre Uoo & J^a- If the averages of are 
bounded, we may assume by compactness that these averages converge; this, however, 
implies that the sequence tt„ itself converges to the element of the ttoo-fibre with limiting 
average. Since this limit is clearly a critical point of F we have the contradiction in this 
case and may assume from now on that the averages of ?7,„ tend nionotonically to +00. 

Let M be such that x > M implies f'{x) > 0. Notice in particular that for all u ^ Si 
there is a t G such that f'{u{t)) = and therefore u{t) < M. For each n, let t„, G S"*^ 
be such that Un{tn) < M. From the compactness of S^, we may assume that the sequence 
tn converges to, say, too- Let Unm ("^ < n) be the element in the ttn-fibre with average 
Um', clearly, Unm < Un- Define similarly Uoom'- the sequence Unm (for fixed m) tends to 
Uoom- Since Unm{tn) < M for all n, Wooml^oo) < M for all m, in contradiction with Lemma 
1.3. □ 

The set S2 of cusps of F is a smooth contractible submanifold of codimension 2 of B^. 
From Lemma 3.1, there is a diffeomorphism of to itself taking the sets Si and ^2 to 

51 — Si and S2, respectively and it suffices to show that ^'2 is contractible. We may thus 
use Corollary 3.3: we check that F is 2-regular, i.e., that is a regular value for both Si 
and (Ei,E2)- We have already seen that DT,i = -D(Si) is never zero in the critical set. 
Also, 

D{ti,±2){u)-v= (^j f'{u{t))v{t)dt,J r{u{t))v{t)dt 

and we are left with showing the linear independence of the functions f"{u{t)) and f"'{u{t)) 
for u satisfying J f'{u{t))dt = J f"{u{t))dt = 0. The first function is non-zero but of 
average zero and the second is strictly positive. □ 

In particular, from Corollary 3.3, there is a homeomorphism of B^ to itself taking 
and S2 to nested subspaces of codimensions 1 and 2. This homeomorphism, however, does 
not respect fibres: we now show how to do better. It is convenient from this point on to 
work in adapted coordinates, i.e., to consider F : B^ B^, its critical set Si = \I/(5'i) 
and set of cusps S2 = \E'(<S'2). Recall that \E' : — > takes fibres to vertical lines and 
that F is in adapted coordinates: in other words, vertical lines are fibres for F. Similarly, 
5° = U 9^6 U He, the disjoint images under of J^a, and J^c- 

In the next claims, we construct a number of auxiliary changes of variable, leading 
eventually to the global normal form for the cusp. In the H case, all constructions are 
smooth. In the C case, however, we make use of homeomorphisms and folds and cusps 
have to be interpreted topologically. Figure 5.2 may be helpful. 

There is a homeomorphism Ti of B^ to itself taking vertical lines to vertical lines and 
3^1, to a subspace C of codimension 1 of B^ . In the H case, this homeomorphism can be 
taken to be a diffeomorphism. From the local normal form for cusps, vertical lines intersect 

52 transversally: the vertical projection is then a natural diffeomorphism between S2 and 
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Figure 5.2 

9^6. In particular, C 5° is a contractible submanifold of codimension 1 for which Lemma 
5.3 apphes: smoothness guarantees the existence of focal tubular neighbourhoods which 
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can be consistently glued because the complement of has two connected components. 
Thus, there is a homeomorphism Tq : taking 5";, to a closed subspace C of 

codimension 1. Define Ti : B^ B^ as the only extension of Tq respecting vertical lines 
and horizontal hyperplanes — Ti is clearly a homeomorphism. □ 

Notice that the conjugation Fi = Ti o F o T^^ is still in adapted coordinates, i.e., 
vertical lines are invariant. Also, the vertical projection of Ti(S2) and Ti(F(S2)) are both 
equal to C; Ti(S2) is the set of cusps (or, in the C case, topological cusps) of Fi and 
Ti(F(S2)) is its image. 

Identify 

s° = ce(r)e(i), 

for r E B^ not in C and now write a typical element of B^ as (Z, x, y) e C x making use 
of the natural projections. Notice that the sign of x determines the type of the (Z, x, 0)- 
fibre: without loss, the cases x < 0, x = and x > are set to correspond to fibres of 
types (a), (b) and (c), respectively. 

There are homeomorphisms T2,T3 of B^ to itself keeping each vertical line invari- 
ant and such that T2(Ti(S2)) = T3(Ti(F(S2))) = C. In the H case, these maps are 
diffeomorphisms. Set T2(Z, x, y) = {Z,x,y — y') where y' is the only real number such 
that TJ~^(Z, 0,y') G S2. Similarly, set T3(Z, x,?/) = {Z,x,y — y") where y" satisfies 
Tr^(Z,0,y")eF(S2). □ 

The composition F2 = TsoFioT^''^ is almost in normal form: for each Z, F2 restricted 
to the Z-plane is a global 2-dimensional cusp. We are now ready to apply Lemma 5.4 to 
get two further changes of variable T4 and T5 such that F3 = T5 o F2 o TJ^ is the desired 
normal form 

F3(Z, X, y) = (Z, X, y^ + xy). 

■ 

Sketch of proof of Lemma 5.4: Item (a) is straightforward. As to item (b), we begin 
by invoking Whitney's construction ([W]) which obtains, for any given Z, diff'eomorphisms 
vz,6 and vz,7 of neighbourhoods of the origin taking the restriction Gz{x, y) = {x, gz{x, y)) 
to the normal form k{x, y) = (x, y^ + xy), i.e., k = vzj o Gz ° '^~z\- The diffeomorphisms 
vz,Q and vzj so constructed are of the form (x, y) 1— > (x, y') and preserve orientation. 

Actually (and here we omit the verification), the construction allows Vzfii and 
the size of the neighbourhoods to be chosen smoothly as functions of the parameter Z. 
More exactly, we have diffeomorphisms g ^ind 7 defined on tubular neighbourhoods 
of Z X (0,0) taking G to the normal form K{Z,x,y) = {Z,x,y^ + xy) near Z x (0,0), 
i.e., K = Ttj o G o g whenever the right hand side is defined. Now extend 5 and 
Ttj to diffeomorphisms Tq and T7 from 2 x to itself of the form (Z, x, y) ^ (Z, x, y') 
which coincide with the identity outside a tubular neighbourhood of Zx (0, 0) — notice that 
this extension is just a one-dimensional problem, parametrized by (Z, x). The composition 
Gi = T7 o G o Tg ^ satisfies all the original conditions in Lemma 5.4 and coincides with 
the normal form K in a tubular neighbourhood of Z x (0, 0). The hard part of bringing 
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the cusps into normal form being done, we give explicit instructions to take Gi to normal 
form. 

Let Sgi and Sk be the critical sets of Gi and K. Their images Gi{Sgi) and K{Sk) 
have two points in each vertical line (Z, x,-), x > 0. Construct Tg by juxtaposing 1- 
dimensional maps to be a diffeomorphism of Z x (0, 0) satisfying 

• T9(Z,a;,y) = {Z,x,y'), 

■ for each Z there is a positive xz such that Tq{Z, x, y) = {Z, x, y) for x < xz, 

■ T^{Gi{Sg,))=K{Sk). 

Let ^2 = ^906*1; notice that G2{Sg2) = K{Sk), where S'g^ ^^e critical set of 6*2- 
There is now a unique diffeomorphism Tg satisfying G2 o Tg ^ = K. Indeed, if x < 0, let 
TsiZ, X, y) be the only point (Z, x, y') for which G2{Z^ x, y) = K{Z, x, y'). For positive x, 
given (Z, x, y) there are at most three points satisfying G2{Z, x, y) = K{Z, x, y'): we define 
Ts{Z, X, y) to be the only point {Z, x, y') which is in the same position with respect to the 
two critical points of K in the vertical line (Z, {Z, X, y) is with respect to the two 

critical points of G2 in the same vertical line. Clearly, Tg is a homeomorphism: we check 
smoothness of Tg and its inverse. At regular points, this is the inverse function theorem. 
At fold points, one may use the square root trick in the proof of Theorem 4.2, but we omit 
the details. Finally, smoothness at cusps is guaranteed from the simple fact that Tg turns 
out to be the identity near Z x (0, 0). ■ 

As a corollary, we obtain a global cusp form for the Cafagna-Donati operator ([CD], 
[CT]): 

Corollary 5.5: Let f{x) ~ ax + bx^ + cx'^'^^^ where k is a positive integer, a > 0, 
+ 6^ > and c < 0. Then the operator F : —>■ is a smooth global cusp and 
F : —>■ is a topological global cusp. 

6. A numerical counter-example 

It is of course tempting to speculate about the possible consequences of D^/ > 0: does 
this condition at least guarantee that points have at most four pre-images? In [L], Lins 
Neto shows that, if f{x,t) is a polynomial of degree four in x with coefficients depending 
on t and positive highest degree coefficient, then the number of solutions may be arbitrarily 
large. 

In this section, we obtain a polynomial / of degree 4 and a smooth periodic Ub such 
that Ub is a Morin singularity of order 4 (a buttertiy). From Morin's normal form, some 
points V near F{ub) have five regular pre-images close to u^. Since the degree of F is zero, 
there is yet a sixth pre-image and we thus obtain a smooth periodic function v for which 
the equation 

u' + f{u)=v{t),u{Q) = u{l), 

has six solutions. 
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We briefly describe the numerical procedure used in the identification of Ub. Without 
loss of generality, f{x) = — bx^ + cx and we try to find a butterfly ui of the form 

Ub{t) = ao + ai cos(t) + a2 cos(2t) + 62 sin(2t) + h ^4 sin(2t). 

We now write the four scalar equations T,i{ui,) = 0, i = 1,...,4, in terms of the ten 
parameters 6, c, ao, . . . , 64 and search for a zero with a Newton-like method with pseudo- 
inverses [AG] . Actually, for appropriate b and c, four extra parameters should be enough 
to locate a butterfly, but the numerical analysis becomes more robust with additional 
parameters. One example is 6 = 4, c = —0.3, oq = —0.01173378, ai = —0.8836063, 
02 = 0.2428734, 62 = -0.6855379, 03 = 0.4465347, 63 = 0.1853376, 04 = -0.01881213 and 
64 = 0.2105862. The Newton method itself checks for the surjectivity of (the restriction 
of) the derivative of (Si, . . . , E4) and the program also verifles that T,5{ui)) ^ 0. 




Figure 6.1 



Again by Newton's method, we try to solve 

(El, ... , E4)(wi) = (-0.0000005, 0, 0.00008, 0), 

where the non-zero constants on the right hand side were adjusted somewhat empirically — 
in a nutshell, we are trying to perturb the polynomial to get five distinct real roots, 
which can be accomplished by adding small multiples of x^ and x. The parameters for ui 
are 

ao = -0.011367708203969, ai = -0.883600656945802, 
02 = 0.243308077825844, as = 0.446085678376277, 04 = -0.018458472190807, 
62 = -0.685621717642052, 63 = 0.185481811055651, 64 = 0.210509692732880. 
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In Figure 6.1 we plot Pv{x) — x, where v = u[ + f{ui) = F{ui), so that roots of this 
auxiUary function correspond to periodic solutions of u' + f{u) = v. This graph was 
obtained by solving the differential equation with a Runge-Kutta method of order 4 for 
initial conditions ranging from —0.4 to 0.4. The vertical scale is stretched by a factor 
of 2.5 • 10^ and the time step for the method had to be taken as 2 • 10~^. Notice the 
clustering of the first five roots, stemming from the butterfly: actually it is this clustering 
and the quintic behaviour of the butterfly which account for the need of a huge vertical 
stretching factor. Another consequence of the quintic behaviour is the great sensitivity 
of the coefficients: for instance, a change of 10~^ in gq destroys four of the six solutions. 
Still, this final direct check is far easier (and more reliable) than the process of obtaining 
the coefiicients for the example. 
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